Open Access 



Research 



bmj Individual-level, network-level and 

0[D6n city-level factors associated with 
HIV prevalence among people 
who inject drugs in eight Russian 
cities: a cross-sectional study 



Ksenia Eritsyan, 1 ' 2 Robert Heimer, 3 Russell Barbour, 3 Veronika Odinokova, 1 
Edward White 3 Maia M Rusakova, 1 Tatiana T Smolskaya, 4 Olga S Levina 1 



To cite: Eritsyan K, Heimer R, 
Barbour R, etal. Individual- 
level, network-level and 
city-level factors associated 
with HIV prevalence among 
people who inject drugs in 
eight Russian 
cities: a cross-sectional 
study. BMJ Open 201 3;3: 
e002645. doi:10.1136/ 
bmjopen-201 3-002645 

► Prepublication history for 
this paper are available 
online. To view these files 
please visit the journal online 
(http://dx.doi.org/10.1136/ 
bmjopen-201 3-002645). 



Received 28 January 2013 
Revised 3 April 2013 
Accepted 5 April 2013 



This final article is available 
for use under the terms of 
the Creative Commons 
Attribution Non-Commercial 
2.0 Licence; see 
http://bmjopen.bmj.com 



For numbered affiliations see 
end of article. 



Correspondence to 

Ksenia Eritsyan; 

ksenia.eritsyan@gmail.com, 

ksenia@ngostellit.ru 



ABSTRACT 

Objectives: To ascertain HIV prevalence among 

people who inject drug (injection drug users (IDUs)) in 

the Russian Federation and identify explanations for the 

disparity in different cities. 

Design: Cross-sectional survey with serological 

testing for HIV and hepatitis C virus prevalent 

infections. 

Setting: 8 Russian cities— Irkutsk, Omsk, 
Chelyabinsk, Yekaterinburg, Naberezhnye Chelny, 
Voronezh, Orel and St Petersburg. 
Participants: In 2007-2009 active IDUs were 
recruited by respondent-driven sampling with a target 
sample size of 300 or more in each city. 
Main outcome measures: Participants were 
administered a questionnaire covering 
sociodemographics, injection risk and protective 
behaviours, sexual behaviours, HIV knowledge, 
experiences with drug treatment and harm reduction 
programmes and social networks. Participants were 
tested for HIV and hepatitis C by enzyme 
immunoassay. Data were analysed to identify 
individual-level, network-level and city-level 
characteristics significantly associated with HIV 
prevalence. Factors significant at p<0.1 were entered 
into a hierarchical regression model to control for 
multicollinearity. 

Results: A total of 2596 active IDUs were recruited, 
interviewed and tested for HIV and hepatitis C virus 
infection. HIV prevalence ranged from 3% (in 
Voronezh) to 64% (in Yekaterinburg). Although 
individual-level and network-level variables explain 
some of the difference in prevalence across the eight 
cities, the over-riding variable that seems to account 
for most of the variance is the emergence of 
commercial, as opposed to homemade, heroin as the 
predominant form of opioid injected. 
Conclusions: The expansion of commercial heroin 
markets to many Russian cities may have served as a 
trigger for an expanding HIV epidemic among IDUs in 
that country. 



ARTICLE SUMMARY 



Article focus 

■ To test the hypotheses that differences in HIV 
prevalence among injection drug users (IDUs) 
could be explained by different level variables 
such as individual injection practices, the type of 
drug(s) injected, average duration of drug use 
structure of IDU populations, interaction within 
social networks of IDUs or a combination of 
these factors. 



Key messages 

■ The difference in HIV prevalence is best 
explained by city-level difference in the type of 
drug injected. 

■ HIV-prevention programmes targeting injectors 
that increase knowledge and reduce injection risk 
may reduce transmission, even in cities that have 
already experienced large epidemic. 

Strengths and limitations of this study 

■ The largest side-by-side comparison of HIV and 
hepatitis C virus prevalence among IDUs linked 
to individual level behavioural and higher order 
social network and city level data yet conducted 
in Russia. 

■ Change in the structure of the heroin market 
emerged as a possible explanation for the differ- 
ences in HIV prevalence detected in the eight 
cities. 

■ Non-random selection of the eight cities, non- 
probability respondent-driven sampling method. 

■ Cross-sectional design prevents determination of 
temporality needed to demonstrate causative link 
between changing heroin market and increasing 
HIV prevalence. 



INTRODUCTION 

The Russian Federation has experienced one 
of the fastest growing HIV/ AIDS epidemics 
observed anywhere in the world since 
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significant spread began in 1996. 1 2 The number of offi- 
cially registered people with HIV increased over 250-fold 
between 1996 and 2011, from 2600 to more than 
650 000 officially diagnosed cases. 3 Most recent estimates 
place the number of infections countrywide at around 
1 million (840 000-1 200 000), with an estimated HIV 
prevalence among Russian adults in 2009 around 1% 
(0.9-1.2%). 2 

Since 1996, the HIV epidemic in Russia has been 
largely driven by injection drug use. There are almost 
400 000 officially registered injection drug users (IDUs) 
in Russia 5 ; the estimated number is several times higher: 
1.5-2.2 million. 6-8 IDUs comprise nearly 80% of the 
registered HIV cases. 4 Surveillance data indicate great 
regional differences in HIV prevalence among IDU 
population; the rates of HIV-prevalence among IDUs in 
governmental treatment and rehabilitation programmes 
range between 0% and 48%. 5 These data are supported 
by research that found HIV prevalence among IDUs 
ranging from 3% to 56%. 9-11 Factors that may give rise 
to the observed variation in HIV prevalence include 
individual-level parameters (eg, hygienic or unhygienic 
injection and sex practices), social network parameters 
(eg, factors that may influence the proportions of IDUs 
in network who are exposed) and city-level parameters 
(eg, time of entry of the virus into the local IDUs popu- 
lation, size of IDU population, type of drug and charac- 
teristics of drug markets) . 

By 2012, the official case notification surveillance data 
at the outset of the study from the oblasts (provinces) in 
which these eight cities are located revealed a wide 
range in HIV prevalence for the province as a whole — 
from less than 40 cases/ 100 000 in Voronezh oblast 
(38.6) to more than 1000 cases/100 000 in Sverdlovsk 
oblast (Yekaterinburg, 1063.5) and Irkutsk oblast 
(1349.3), with intermediate rates in the other regions. 3 
In addition to HIV, hepatitis C virus (HCV) is endemic 
among IDUs. 12 " 14 However, data on the prevalence of 
HCV in Russian IDUs are insufficient. In the few 
Russian locales where serological data on HIV and HCV 
are available, HCV prevalence exceeds that for HIV 11 15 
This is not only true in Russia but also in most of the 
world. 16 

Major differences among IDUs within Russia are the 
type and formulation of the drugs they inject. The two 
dominant drugs are heroin and methamphetamine-type 
stimulants (MTS) . Both have a history of home manufac- 
ture — heroin from poppies or opium gum and MTS 
from ephedrine, pseudoephedrine or phenylpropanola- 
mine — and are frequently prepared, stored and sold to 
IDUs as liquids. 17-20 Homemade MTS involves reduction 
to methamphetamine or amphetamine, or oxidation to 
methcathinone or cathinone. Since our ethnographic 
studies of Russian IDUs who manufacture or use MTS 
suggested that they often confused oxidation and reduc- 
tion and could not reliably report what MTS they either 
made or injected, 18 21 we group all these drugs into a 
single MTS category. Homemade heroin involves the 



extraction of the morphine in alkaline, organic extrac- 
tion, acid back extraction and conversion of the mor- 
phine to heroin (diacetylmorphine). However, in many 
parts of Russia, homemade heroin has been replaced by 
commercial heroin trafficked directly from Central 
Asia. 7 22 This transition took place in the late 1990s 
jointly with a rapid increase of HIV epidemic in Russia. 23 
Our interest in the variation in HIV prevalence among 
IDUs was engaged by the opportunity to conduct 
behavioural and serological data collection and analysis 
in eight diverse Russian cities — Irkutsk, Omsk, 
Yekaterinburg, Chelyabinsk, Naberezhnye Chelny, 
Voronezh, Orel and St Petersburg using the same sam- 
pling approach — respondent-driven sampling (RDS) — 
and virtually identical data collection instruments in all 
eight cities. The work has yielded the largest side-by-side 
comparison of HIV and HCV prevalence linked to 
individual-level behavioural and higher order social 
network and city level data yet conducted in Russia. 
Observations from prior studies 19 24-26 led us to formu- 
late the hypotheses that differences in HIV prevalence 
will be explained by different level variables such as: 
individual injection practices, the type of drug(s) 
injected, average duration of drug use, percentage of 
men and women in particular IDU populations, inter- 
action within social networks of IDUs or a combination 
of these factors. The resulting analysis permits an explor- 
ation of the associations of HIV infection with 
individual-level, network-level and city-level parameters. 
Owing to the factors that influence the likelihood of 
infection acting on different levels, it is necessary to 
develop a hierarchical approach in forging analytical 
strategies. 

METHODS 

Recruitment procedure 

Recruitment using RDS of IDUs was conducted in eight 
Russian cities — Voronezh, Orel, Naberezhnye Chelny, 
Yekaterinburg, Chelyabinsk, Omsk and Irkutsk — with a 
target of 300 IDUs in each city and St Petersburg with 
target sample of 400 IDUs (figure 1). The sample size 
was set by the study hinders and based on many inter- 
national studies they were supporting at the same time. 
For inclusion in the study, participants had to be at least 
15 years old and to report having injected drugs in the 
last month. In accordance with RDS methodology, the 
initial participants of the study, known as 'seeds', were 
members of the target IDU population in their respect- 
ive cities, and were selected based on their willingness to 
recruit their peers into the study. 27 28 Sample sizes were 
set to achieve samples consistent with previous work that 
held such sample sizes reached sample equilibrium and 
reduced the biases introduced by non-random selection 
of seeds. 29-31 Recruitment was initiated with 6-7 seeds in 
each city referred from local organisations, governmen- 
tal and/or non-governmental, conducting an outreach 
or the other kind of work with the IDU population. 
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Figure 1 Recruitment using respondent-driven sampling of injection drug users was conducted in eight Russian cities- 
Voronezh, Orel, Naberezhnye Chelny, Yekaterinburg, Chelyabinsk, Omsk and Irkutsk. 



Non-productive seeds were replaced such that in sum 
the sample was collected relying on 6-24 productive and 
0-75 unproductive seeds in each city. A total of 2582 
IDUs were recruited between November 2007 and 
October 2009. 

Each person recruited and enrolled in the study 
received three recruitment coupons to recruit his/her 
peers in the study. Coupon distribution ceased when the 
sample size was reached after which participants were 
warned that the study would be ending within a few 
days. Coupons consisted of three parts that could be 
easily separated. All parts had the same unique number 
printed on them. One part of the coupon served as the 
'referral coupon', which the recruiter used to recruit his 
peers into the study and the other part served as the 
'payment coupon' and was kept by the recruiter to claim 
an incentive for having recruited his peers into the 
study. The third one has to be used by the respondent 
wanting to know the results of his/her blood testing. 
Each enrolling participant was given a unique consecu- 
tive ID number in addition to the coupon number. 

Participants were given a gift package consisting of 
food (chocolate, a can of coffee or tinned goods) and 
personal hygiene products (shaving items and 
shampoo). The gifts were equivalent to US$14 for par- 
ticipation and US$7 for each recruit referred to and 
enrolled into the study. Recognition of potential dupli- 
cate participants depended upon study staff who were 
outreach workers familiar with the target population. 

The final sample sizes of participants for whom both 
behavioural and serological data were collected in each 
city were 279 in Irkutsk, 350 in Omsk, 306 in 



Chelyabinsk, 300 in Yekaterinburg, 341 in Naberezhnye 
Chelny, 309 in Voronezh, 300 in Orel and 411 in 
St Petersburg. Sample accrual time ranged between 5.5 
and 22 weeks. 

Collection of demographic, behavioural and 
serological data 

The data collection was preceded by a written informed 
consent. To ensure the participants' anonymity, they 
were allowed to put a false signature on the consent 
paper. Data collection consisted of face-to-face interviews 
and blood specimens. The data collection team at each 
site consisted of one nurse, one or two interviewers and 
one person serving as coupon manager, all supervised 
by a study coordinator. The coupon manager recorded 
data for each participant on a form that included each 
individual ID, the coupon number brought to the inter- 
view and the coupon numbers each new participant was 
given for distribution. Interviewers were chosen based 
on their personal experience working with IDU popula- 
tions (social workers, psychologists or epidemiologists 
from the local NGO or AIDS Center). 

In Naberezhnye Chelny, Voronezh, Yekaterinburg and 
Omsk data collection occurred at a single site based in 
Regional AIDS Center and in Chelyabinsk and Orel data 
collection occurred at a single site operating from a van 
parked at one fixed location. In Irkutsk, there were mul- 
tiple study sites: one site was based in Regional AIDS 
Center and two or three others were located outside the 
city centre and operated special hours at weekdays from 
mobile vans at fixed locations. In St Petersburg, five sites 
operated at regional narcology (drug treatment) clinics. 
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After the interview and pretest counselling, blood 
specimens were collected by venipuncture (3-5 ml of 
blood) using closed vacuum systems into tubes marked 
with the participant's unique ID (a combination of 
coupon and consecutive participation number). HIV 
testing was performed using an enzyme-linked immuno- 
assay with confirmatory western blot at the laboratories 
of Regional AIDS Center in each participating city. 
Participants were instructed to return in 7-14 days for 
their test results. During the post-test visit, the 
HIV-positive participants were counselled and referred 
for long-term care to physicians working at Regional 
AIDS Centers. 

Questionnaire topics included network sizes informa- 
tion, sociodemographic indicators, sexual and drug 
injection risk behaviours during the 30 days prior to 
interview, knowledge of HIV and STIs and previous HIV 
testing. The questionnaire was based on FHI guidelines 
for conducting behavioural surveillance. 32 Survey dur- 
ation was approximately 1 h. 

Data analysis 

Data from the eight cities were compiled into a single 
database. Descriptive statistics, bivariate associations 
between HIV prevalence and other variables and hier- 
archical models were produced using SPSS and R soft- 
ware. Initially, three models to identify the factors 
associated with HIV prevalence were generated. Because 
missing data represented <4% (106 of 2596) of indivi- 
duals, they were excluded from the models. In the first 
model, only individual-level variables with bivariate asso- 
ciations with p<0.1 were entered into a backward step- 
wise process using the 'step' command in R, which is 
based on the best Akaike Information Criteria. In the 
resulting multivariate model, all remaining variables 
were considered significant if p<0.05. A second model 
added to the first model three measures of network 
characteristics — participants' self-reported network size 
(number of IDUs participant met personally during the 
last 6 months), the mean number of participants in the 
RDS chain in the city and the frequency of meeting 
one's recruiter. This second model contains all those 
variables that remained significant at p<0.05. A third 
model added into the second model three city level vari- 
ables. The first two variables — the dominant form of 
heroin (homemade vs commercial) and high or low pro- 
portions of individuals injecting MTSs — are dichotom- 
ous; the third — the average length of RDS chains in 
each city — is a continuous variable. Again, the third 
model contains all those variables that remained signifi- 
cant at p<0.05. 

The multicollinearity of the city-level variable of pre- 
dominant commercial heroin and all other variables 
significant in the third model was tested using bivariate 
analysis. A binomial distribution within generalised 
linear models (GLM) for both univariate and multivari- 
ate analyses was implemented using R software (http:// 
r-forge.r-proj ect.org). Hierarchical modelling applied 



the <hglm> command within the hglm package in R, 
specifying commercial availability of heroin by city as a 
random variable, with demographic, network and 
individual injection risk behaviours that were found stat- 
istically significant in the aforementioned GLM specified 
as fixed variables. 



RESULTS 

A comparison of the demographic data from across the 
eight cities revealed strong similarities for seven of the 
eight cities (table 1). In these seven cities, the percent- 
age of male IDUs ranged from 69% to 76% and the 
mean and median ages were in the range 26-31 and 
26-30, respectively. The exception was Chelyabinsk, 
where 89% of IDUs were men and the mean and 
median ages were 24.7 and 24, respectively. In addition, 
the duration of injection among IDUs in Chelyabinsk 
was 6 years, less than in any city other than Orel. The 
majority of participants had completed secondary educa- 
tion through the Russian equivalent of the 10th grade, 
and 7-27% in each city had at least some postsecondary 
education. 

HIV prevalence, the primary dependent variable, 
ranged from 3% in Voronezh to 64% in Yekaterinburg. 
HCV prevalence ranged from 49% in Naberezhnye 
Chelny to 90% in Yekaterinburg and St Petersburg. 

In terms of individual-level injection risk behaviours, the 
majority of participants in each city reported never having 
been injected with a used syringe during the 30 days prior 
to interview (range 61-86%). Many, however, during this 
same period reported sharing other injection parapherna- 
lia (range 26-91%) and engaging in other practices with 
the potential to transmit HIV or HCV (table 1). 
Engagement with services for IDUs was sporadic. Most had 
no recent contact with harm-reduction programmes 
(range 65-95%) and only one-quarter were in a drug treat- 
ment programme in the preceding year (range 7-45%). 

The drug injected most frequently was one of several 
forms of heroin — either commercially trafficked or 
home manufactured (from opium, poppy straw or 
poppy seeds) . The data in table 1 reveal that most IDUs 
in any given city inject with either the homemade or 
commercial variety, but not both. Prevalence of injection 
of commercial heroin was 95-100% in six of the eight 
cities — Irkutsk, Omsk, Chelyabinsk, Yekaterinburg, 
Naberezhnye Chelny and St Petersburg — but uncom- 
mon in the other two cities. Injection of MTS was less 
common, with prevalences ranging from 30% in 
Yekaterinburg to 2% in Voronezh and 0% in Irkutsk. 

We developed stepwise models to identify factors asso- 
ciated with HIV prevalence across the eight cities. In the 
first model, which included only the individual-level vari- 
ables from table 1, we found that 7 of the 14 were asso- 
ciated in bivariate analysis (table 2), but when these 
seven were included in a multivariate model only four 
remained significant. Injecting with someone known to 
be HIV-positive and sharing rinse water were positively 
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179 (58%) 


1 1 1 (37%) 


86 (21%) 


Injected methamphetamine-like stimulants 


0 


14 (4%) 


30 (10%) 


89 (30%) 


62 (18%) 


7 (2%) 


79 (26%) 


30 (7%) 



Table 1 Continued 



Naberezhnye 





Irkutsk 


Omsk 


Chelyabinsk 


Yekaterinburg 


Chelny 


Voronezh 


Orel 


St Petersburg 


Individual level protective factors 


















Injected with sterile syringe 


















Always in last 30 days 


177 (63%) 


215 (61%) 


192 (63%) 


241 (80%) 


229 (67%) 


193 (62%) 


257 (86%) 


280 (68%) 


Harm reduction services 


















None 


234 (84%) 


294 (84%) 


198 (65%) 


225 (75%) 


276 (81%) 


271 (88%) 


284 (95%) 


361 (88%) 


<Once/month 


13 (5%) 


27 (8%) 


17 (6%) 


17(6%) 


10 (3%) 


7 (2%) 


7 (2%) 


32 (8%) 


>Once/month 


20 (7%) 


28 (8%) 


90 (29%) 


59 (20%) 


52 (15%) 


31 (10%) 


7 (2%) 


18 (4%) 


Missing data 


12 (4%) 


1 (0%) 


1 (0%) 


0 


3(1%) 


0 


2 (1%) 


0 


Governmental drug treatment 


















Ever 


162 (58%) 


60 (17%) 


89 (29%) 


45(15%) 


204 (60%) 


105 (34%) 


63 (21%) 


152 (37%) 


Last year 


126 (45%) 


26 (7%) 


49 (16%) 


24 (8%) 


146 (43%) 


93 (30%) 


27 (9%) 


1 1 1 (27%) 


Perfect HIV knowledge 


145 (52%) 


123 (35%) 


199 (65%) 


108 (36%) 


204 (60%) 


198 (64%) 


66 (22%) 


239 (52.8%) 


Network level variable 


















RDS network size, past 6 months 


















Mean 


35.8 


39.9 


10.8 


49.7 


16.6 


21.9 


17.3 


10.4 


Median 


15 


20 


9 


22 


10 


15 


10 


5 


Range 


1-500 


1-1000 


2-50 


1-1000 


1-500 


1-200 


1-150 


1-500 


City level variables 


















Commercial Heroin Dominant (% of IDUs 


Yes (98%) 


Yes 


Yes (98%) 


Yes (99%) 


Yes (100%) 


No(1%) 


No (31%) 


Yes (95%) 


reporting 




(99%) 














any use, past 30 days) 


















Homemade heroin Dominant (% of IDUs 


No (6.5%) 


No (7.4%) 


No (6.5%) 


No (19.7%) 


No (2.6%) 


Yes 


Yes 


No (0.5%) 


reporting any use, past 30 days) 












(96.8%) 


(79.0%) 




>1 in 6 IDUs reported Recent MTS use 


No (0%) 


No (4%) 


No (10%) 


Yes (30%) 


Yes (18%) 


No (2%) 


Yes (26%) 


No (7%) 


RDS chain length 


















Mean 


33.0 


100.6 


112.0 


55.7 


60.8 


83.8 


98.9 


39.9 


Median 


33 


86 


169 


55 


48 


83 


156 


35 


Range 


1-65 


1-164 


5-169 


1-96 


23-90 


1-120 


1-156 


1-85 



CD 



1 



HCV, hepatitis C virus; IDU, injection drug user; RDS, respondent-driven sampling. 
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Table 2 Association of individual-level demographic, risk and protective factors with HIV prevalence in eight Russian cities 


Variable 


Crude OR (±95% CI) 


Adjusted OR (±95% CI) 


bex 


1 .09 (O.yO to 1 .31) 




Age 


1 .01 (0.99 to 1 .02) 




Educational level 


U./o (U.oo TO 1 .1 /) 




Duration of injection 


1 .UZ (1 .01 TO 1 .04) 


a r\A ir\ r\r\ t. n a r\0\ 

i .01 (o.yy to i .03) 


Injected with used syringe, last 30 days 


^ AG. /H OH H ~?G\ 

1 .4b (1 to 1 . /b) 


U.ob (U.bb TO 


Injected with known HIV+, last 30 days 


A C\r^ IA ~7 A 4- « O OH\ 

1 .96 (1 .74 to 2.21) 


1 .86 (1 .64 to 2.1 1) 


Shared paraphernalia, last 30 days 


1 .58 (1 .33 to 1 .88) 


A A a / r\ nt~\ X A o o \ 

1.11 (0.89 to 1 .38) 


Shared rinse water, last 30 days 


o on /h no o ~7A\ 
d.o\) (1 .yo TO /4) 


O OQ /H "7Q O QC\*** 


Injected with prefilled syringes, last 30 days 


1.06 (0.88 to 1.28) 




Injected methamphetamine-like drugs, last 30 days 


1.10 (0.88 to 1.39) 




Always injected with sterile syringe, last 30 days 


0.99 (0.82 to 1.19) 




Used harm reduction services, last 12 months 


0.77 (0.61 to 0.98) 


0.74 (0.56 to 0.97)* 


Ever in treatment for addiction 


0.92 (0.76 to 1.10) 




Perfect score on HIV knowledge test 


0.73 (0.62 to 0.87) 


0.78 (0.65 to 0.95)* 


***p<0.001, **p<0.01, *p<0.05. 



associated, whereas the recent use of a syringe exchange 
programme and perfect score on the HIV knowledge 
test was protective. None of the sexual risk behaviours 
considered in the survey instrument (unprotected 
sexual intercourse, number of recent or HIV-positive 
partners and engagement in commercial sex) was signifi- 
cantly associated with HIV prevalence. 

In the second model, we determined which of the three 
network-level variables (RDS network size, mean chain 
length and frequency of meeting the recruiter) were cor- 
related with HIV prevalence and added these to the four 
individual-level variables (table 3). Although both larger 
self-reported network size and shorter length of an indivi- 
dual's referral chain were significantly associated with 
prevalence in bivariate analysis, the multivariate model 
excluded the length of the referral chain and did not elim- 
inate any of the four significant individual-level variables. 

In the third model, we determined which of the three 
city-level variables were correlated with HIV prevalence 
and added these to the five previously identified variables 
(table 4). All were found to be significantly associated with 
HIV prevalence in bivariate analysis, but only two were in 
the multivariate analysis. The predominance of commer- 
cial heroin was positively associated with prevalence, 
whereas larger mean recruitment chain length was pro- 
tective. One of the individual-level variables, recent use of 



syringe exchange, was no longer protective in this model. 
Of the six factors associated with HIV prevalence, the vari- 
able with the strongest predictive value was the city-level 
variable commercial heroin predominance. We explored 
the possibility of autocorrelation among these variables as 
well as variables significantly associated with HIV preva- 
lence in all previous models (table 5). We found that of 
the eight variables, only two — sharing cookers and higher 
levels of MTS injecting — were not correlated and all the 
variables that remained significant in the third model were 
collinear with commercial heroin predominance. 

Given the high degree of collinearity, we constructed 
one final model using hierarchical modelling controlling 
for commercial heroin predominance as a random effect 
(table 6). The results of this model are consistent with the 
third model presented in table 4 despite the high collin- 
earity of the significant variable. It is interesting to note 
that in the hierarchical model, the intercept is not statistic- 
ally significant, suggesting that incorporating the variable 
of commercial heroin predominance as a random variable 
accounts for much of the variation in the data. 

DISCUSSION 

This study, using bio-behavioural surveillance data on HIV 
prevalence among people who inject drugs from eight 



Table 3 Association of individual-level and network-level 


factors with HIV prevalence in eight Russian cities 


Variable 


Crude OR (95% CI) 


Adjusted OR (95% CI) 


RDS network size, last 6 months 


1.003 (1.001 to 1.004) 


1.003 (1.001 to 1.005)** 


Length of referral chain 


0.998 (0.996 to 0.999) 


0.999 (0.997 to 1.0006) 


Frequency of meeting referrer 


1.09 (0.91 to 1.32) 




Injected with known HIV+, last 30 days 




1.86 (1.64 to 2.1)*** 


Shared rinse water, last 30 days 




2.26 (1 .86 to 2.74)*** 


Used harm reduction services, last 12 months 




0.75 (0.57 to 0.98)* 


Perfect score on HIV knowledge test 




0.76 (0.63 to 0.92)** 


*p<0.05, **p<0.01 , ***p<0.001 . 
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Table 4 Association of individual-level, network-level and city-level factors with HIV prevalence in eight Russian cities 


Variable 


Crude OR (95% CI) 


Adjusted OR (95% CI) 


Mean city recruitment chain length 
L/urrirrieruicii neroin preuoiuinani in uiiy 
Amphetamine-type drugs more common 
RDS network size, last 6 months 
Injected with known HIV+, last 30 days 
Shared rinse water, last 30 days 
Used harm reduction services, last 12 months 
Perfect score on HIV knowledge test 


0.99 (0.98 to 0.99) 

D. 1 1 ^4.00 IO O.c.H) 

1.29 (1.09 to 1.54) 


0.992 (0.989 to 0.995)*** 

\*L.OO IO H-.ODj 

1.17 (0.94 to 1.45) 
1.0025 (1.001 to 1.004)** 
1.58 (1.39 to 1.80)*** 
1.94(1.58 to 2.37)*** 
0.87 (0.66 to 1.16) 
0.79 (0.65 to 0.96)* 


*p<0.05, **p<0.01 , ***p<0.001 . 



widely separated and divergent Russian cities, provides the 
clearest picture to date of the Russian HIV epidemic in 
this at-risk group. The thrust of the analysis, which is 
exploratory in nature, seeks to identify the characteristics 
at the individual, network and city level that might provide 
some explanation for the differences in HIV prevalence 
across the eight cities in the Russian Federation. Three dis- 
tinct features of the epidemic emerge from the analysis. 
First, among the best explanations for the substantial het- 
erogeneity in HIV prevalence, is a strong association 
between the nature of the drug market and HIV preva- 
lence. Second, is the finding that there appears to be 
about a 20% reduction in the likelihood of HIV infection 
among those individuals with greater HIV knowledge and 
a 25% reduction among those who availed themselves of 
syringe exchange programme (although the latter associ- 
ation disappeared when drug markets were considered). 
Third, on the individual-level shared paraphernalia other 
than syringes (in our case — shared rinse water) is an 
important factor of HIV transmission. Each of these will be 
considered in the paragraphs below High prevalence was 
not restricted to any one region of Russia, but was seen in 
cities ranging from Irkutsk in central Siberia (and the 
eastern and southernmost city in the study) to St 
Petersburg (the western and northernmost city in the 
study). There was also no apparent association between 
population size and HIV prevalence (Spearman r s =0.407, 
p>0.05), suggesting that the epidemic was not simply con- 
centrated in larger cities. Instead, HIV prevalence is 



higher in cities in which the form of heroin has made the 
transition from homemade to commercial. The explan- 
ation for this association is unclear, but several non- 
exclusive possibilities exist. 

First, it is possible that the relative ease of using com- 
mercial heroin — no need to spend time manufacturing 
the heroin or to inject large volumes of weak drug solu- 
tions — has led to increases in the number of IDUs and 
hence to changes in the social networks whose interrela- 
tionship and whose injection behaviours may foster the 
spread of blood-borne viruses. The correlation of 
smaller network sizes across different measures of the 
network size may be a consequence of market changes 
in which fewer connections are needed to obtain com- 
mercial heroin compared with that needed to secure a 
reliable supply of homemade opiates. That many of the 
injection risk and network variables are collinear in 
cities where commercial heroin is dominant (as shown 
in table 5) is supportive of but by no means proof of 
such a hypothesis. Factors such as the structure, density 
and stability of those networks could be the key ones 
that lead to differences in HIV-related risks in among 
IDU networks. 34 35 A second possible explanation lies in 
the differences in practices regarding the use of com- 
mercial and homemade heroin that include how these 
drugs are distributed, prepared and injected. Such phe- 
nomena have been described in other studies suggesting 
that the formulation of heroin can play a significant role 
in drug preparation and administration practices linked 



Table 5 Collinearity of city-level commercial heroin predominance and other variables 


Variable name 


Variable level 


Z statistic 


p Value 


Injected with known HIV+, last 30 days 


Individual 


10.16 


<0.00001 


Shared sinse water, last 30 days 


Individual 


8.192 


<0.00001 


Shared cooker 


Individual 


-1.713 


0.0867 


Used harm reduction Services, last 12 months 


Individual 


-3.354 


<0.00001 


Perfect score on HIV Knowledge test 


Individual 


2.461 


0.0138 


Individuals' RDS network size 


Network 


1.993 


0.0463 


Mean city recruitment Chain length 


City 


-16.62 


<0.00001 


Amphetamine-type drugs more common 


City 


1.233 


0.218 



RDS, respondent-driven sampling. 
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Table 6 Hierarchical Generalised Logistic Model of HIV prevalence in eight Russian cities 




Variable 


Estimate 


t value 


d Value 


Adjusted OR (95% Ch 


Mean city Recruitment chain length 


0.0088 


-4.729 


<0.0001 


0.99 (0.987 to 0.994) 


Amphetamine-type drugs more common 


0.084 


0.685 


0.49 




RDS network size, last 6 month 


0.0025 


-2.898 


0.0038 


1.003 (1.001 to 1.004) 


Injected with known HIV+, last 30 days 


0.483 


7.197 


<0.0001 


1.61 (1.41 to 1.84) 


Shared rinse water, last 30 days 


0.679 


6.368 


<0.0001 


1 .96 (1 .58 to 2.42) 


Used harm Reduction services, last 12 months 


0.140 


0.937 


0.35 




Perfect score on HIV knowledge test 


0.235 


-2.231 


0.026 


0.79 (0.64 to 0.97) 


Intercept 


1.204 


-1.933 


0.053 





to the spread of HIV ~ Laboratory data simulating 
Russian drug preparation practices demonstrated that 
HIV is more easily transmitted when using commercial 
heroin. The replication of homemade heroin produc- 
tion, based upon ethnographic observations of the 
process, revealed that not only did HIV fail to survive 
the manufacturing process but also that something 
extracted during the process reduced the viability of 
HIV by 80% when the homemade drug was stored 
inside syringes contaminated with HIV-infected blood. 36 

There also appear to be modest protective effects from 
increased HIV knowledge and attendance at harm reduc- 
tion programmes. The one-fifth reduction in the odds of 
being HIV positive among those who achieved the highest 
score on the HIV knowledge test suggests an important role 
for expanding HIV educational activities targeting IDUs. As 
detailed in many theories of HIV prevention, knowledge is 
necessary, but not necessarily sufficient to reduce risk beha- 
viours and disease transmission. 41-43 Knowledge must be 
complemented by motivation to reduce risk and access to 
prevention tools. On the other hand, according to the data 
presented herein and from other studies, harm reduction 
programmes in Russia have usually attained insufficient 
coverage to be effective in lowering HIV transmission. 44 45 
This may explain why the degree of protective effect is small 
and becomes non-significant once city-level variables are 
considered: although the use of a local SEP may be protect- 
ive in cities where they exist, their existence itself is a city- 
level effect. There was no evidence of protective effect from 
having been in drug abuse treatment, either ever or within 
the past year. This could be considered as an effect of 
Russian narcology policy, which is focused on achieving 
abstinence but not on harm reduction for those who do 
not achieve sustained remission. 46 

Our study shows that sharing non-syringe injection para- 
phernalia is common among Russian IDUs and it signifi- 
cantly correlates which positive HIV status. This means 
that, along with the need to expand the access to clean syr- 
inges that is still insufficient due low NSP coverage and 
IDUs' negative experience of purchasing syringes from the 

to aH rn 

pharmacies, efforts to reduce sharing of parapher- 

nalia, including rinsing water, are urgently needed. 

There are several limitations to this study Two have to 
do with sampling. First, the selection of the eight cities was 



arbitrary — no selection method was developed to ensure 
that the eight cities were representative of the range of 
characteristics of IDUs, drug markets or prevention activ- 
ities found in the Russian Federation. Thus, there is no 
way to generalise these finding to other Russian cities. 
Second, although the same recruitment method was used 
to accrue samples in each city, the recruitment pattern was 
city-specific. This is reflected in the differences in average 
chain length, which ranged from 33 in Irkutsk to 112 in 
Chelyabinsk. There is no sampling frame for random 
selection within the hidden population of IDUs and no 
evidence that the samples of IDUs are representative of 
the underlying population from which the sample was 
drawn; therefore, they must be considered convenience 
samples. A major limitation is that its cross-sectional nature 
makes the study exploratory rather than being definitely 
able to establish causality. The outcome measure of inter- 
est, HIV prevalence, may have occurred at any time in the 
past, whereas the assessment of individual level risk and 
protective behaviours, network attributes and heroin 
market features were ascertained for the recent past only. 
There is a limitation resulting from the data we considered 
in our models. As our modelling was hypothesis-driven, 
variables concerning risky sexual behaviours were not 
included in the models because they failed to predict HIV 
prevalence in any of the eight cities. A final limitation is 
the absence of data on mortality, incarceration and migra- 
tion of IDU populations. Nevertheless, compared with the 
past research, this was conducted in a way that covered the 
largest range of Russian cities yet studied. 

In summary, despite the above limitations, this report, 
which is based on the application of common recruit- 
ment, screening, survey and testing methods, provides 
the best snapshot of the current HIV epidemic among 
IDUs in Russia. The findings reveal that HIV prevalence 
is very high among members of this risk group in several 
cities and that the factor that goes farthest in explaining 
high prevalence is the predominance of a commercial 
heroin market. If a market transition is a pre-requisite 
for increasing prevalence, then subsequent surveillance 
should be conducted in the cities that have experienced 
the transition, but still have low prevalence and in the 
cities that have yet to make this transition. Furthermore, 
more active prevention programmes are urgently 
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needed, including such measures as the increase of 
HIV-related knowledge, expanded syringe and sterile 
equipment access, appropriate addiction treatment, 
increases in antiretroviral therapy specifically for IDUs 
who remain the risk group with highest incidence and 
prevalence and promotion of safer sex and sexual risk 
reduction among IDUs who are at risk of spreading HIV 
to their non-drug-using sex partners. 
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